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The behavior of the charm and bottom structure functions (F^(x,Q 2 ), i=c,b; k=2,L) at small-* 
is considered with respect to the hard-Pomeron and saturation models. Having checked that 
this behavior predicate the heavy flavor reduced cross sections concerning the unshadowed and 
shadowed corrections. We will show that the effective exponents for the unshadowed and saturation 
corrections are independent of x and Q 2 , and also the effective coefficients are dependent to lnQ 2 
compared to Donnachie-Landshoff (DL) and color dipole (CD) models. 


1. Introduction 

In perturbative quantum chromodynamics (PQCD) 
calculations the production of heavy quarks at HERA 
proceeds dominantly via the direct boson- gluon fusion 
(BGF)(Fig.l), where the photon interacts indirectly 
with a gluon in the proton by the exchange of a 
heavy quark pair [1-10]. In the BGF dynamic, the 
charm(bottom) quark is treated as a heavy quark and 
its contribution is realized by fixed- order perturbative 
theory. As to the measurements of HERA [11-15], the 
charm (bottom) contribution to the structure function 
at small a; is a large fraction of the total, since this value 
is approximately 30% (1%) fraction of the total. This 
behavior directly is related to the growth of the gluon 
distribution at small x. We know that the gluons couple 
only through the strong interaction, consequently the 
gluons are not directly probed in DIS. The only way 
to indirect contribution is via the g^qnqH (q h = c, h) 
transition. This involves the computation of the BGF 
process 7 * 3 — >qnqH- This process can be created when 
the squared invariant mass of the hadronic final state 
has the condition that runs as follows W 2 >4m 2 . 

y h 

In the framework of DGLAP (Dokshitzer- Gribov- 
Lipatov-Altarclli-Parisi) [16-18] dynamics, considering 
the heavy-flavour physics is in the framework of variable- 
flavor-number scheme (VFNS) [19,20]. In this scheme, 
the mass logarithms are resummed through all orders 
into a heavy quark density according to the DGLAP 
evolution equations. All logarithmic terms of the heavy 
flavor Wilson coefficients are obtained by factorization of 
this quantity into the massive operator matrix elements. 
The study of heavy flavor production can be done in 
deep inelastic electron-proton scattering, which was 
investigated experimentally at HERA and recently at 
LHC. These data for heavy quark production, have been 
proposed in the framework of the fixed-flavor-number 
scheme (FFNS) (where only light degrees of freedom are 
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considered as active), that it was calculated for F% and 
F$ [11-15], 

The heavy flavor structure functions Fl(x,Q 2 )(k = 2 ,L 
and i = c,b ) are dependent to the parton distribution 
functions. In the small-a: range, further simplification 
is obtained by neglecting the contributions caused by 
incoming light quarks and antiquarks. This is justi¬ 
fied because they vanish at LO and are numerically 
suppressed at NLO for small values of x, while the 
gluon contribution is a matter at this region [6-7]. 
In axial gauges, the leading double logarithms (i.e. 
ln(Q 2 ) ln(l/a:)) are generated by ladder diagrams in 
which the emitted gluons have strongly been ordered 
with respect to the transverse and longitudinal mo¬ 
menta. The sum of theses momenta predicts that the 
gluon distribution increases as x decreases. Clearly this 
increase cannot go on indefinitely. When the density of 
gluons becomes too large, they can no longer be treated 
as free partons [21-22]. At very small x we expect 
annihilation or recombination of gluons to occur and 
these shadowing corrections give rise to nonlinear terms 
in the evolution equation for the gluon distribution 
function. This picture allows us to write the GLR-MQ 
(Gribov, Levin, Ryskin, Mueller and Qiu) equation for 
the gluon distribution function at small-x [23-24]. We 
expect the gluon correlations at small-a; to tame the 
behavior of the gluon distribution function. Therefore 
we observe that the heavy flavor structure functions 
(HFSFs) and also heavy reduced cross sections behaviors 
are tamed the saturation effects. 

An important point in gluon saturation approach is 
the ^-dependent saturation scale Q 2 (x), where it is the 
critical line between non-linear and linear effects and 
it is an intrinsic characteristic of a dense gluon system 
[25-26]. 

In this paper, we also investigate this non-linear behavior 
for the charm and bottom quarks at small-x related to 
the GLR-MQ evolution equation. Then we will apply 
the geometric scaling parameterization in according with 
the critical line Q 2 = Q 2 s {x). We will do this, because 
the geometric scaling of the dipole cross section gives 
the similar scaling of the quantity a s {Q 2 ) xs ^Qp ^ that is 
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dominant in the charm and bottom structure functions. 
The content of our paper is as follows. In the next 
section we give a summary about heavy quark structure 
functions and color dipole model with starting gluon dis¬ 
tribution along the critical line. Then we will study the 
heavy quark structure functions for Q 2 >Q 2 (x) in sec¬ 
tions 3-5, respectively. In Sec.6 we present the behavior 
of the HFSFs exponents at unshadowed and shadowed 
corrections to the gluon density behavior and also in geo¬ 
metrical scaling. Finally we give our conclusions in Sec.7. 


2. A short theoretical input 

In this section we briefly present the theoretical part of 
our analysis. The reader can be referred to the Refs. [27- 
36] for more details. 

The heavy quark contributions F k (x, Q 2 ,mf) to the pro¬ 
ton structure function at small-a; (where only the gluon 
contributions are considerable) are given by this form 


Fl(x,Q 2 ,m 2 ) = 2e 2 - 
xG( 


s(v 2 


2n 

x 

1-z’ 


pl—x 

L-i 

,. 2 \ 


dzci k (i-z,0 


(i) 


where a = l+4£(£=^), G = xg is the gluon distribution 
function and g, 2 denotes the factorization scale. Here 
C g k are the heavy coefficient functions in BGF at LO 
and NLO analysis and in the NLO analysis 

2 4tt 47T/?! Inlnjn 2 / A 2 ) 

’ /3 0 ln( M 2 /A 2 ) Pq ln(jiyA 2 ) [ ’ 

with (3o = 11 — |n/,/3i = 102 — yn.j (n/ is the number 
of active flavors). 

At small-a:, perturbative QCD predicts an increase in the 
gluon distribution tamed by saturation effects. The phys¬ 
ical picture of this process is provided in the rest frame 
of the proton. In the small x limit, the virtual photon 
splits into a qq color dipole followed by the interaction of 
this dipole with the color fields in the proton. The dipole 
cross section has been defined by [37] 

<?dipoie{x, r) = (j 0 (1 - e~ r2Q ' (a:)/4 ), (3) 


in which r is the transverse separation of the qq pair 
and Q s {x) parameterized as Ql(x) = Qq(xo/x) x GeV 2 . 
The important property of the dipole cross section is its 
geometric scaling (GS), which is a well-known property 
of DIS for small-a: values [38-41]. Therefore the proton 
cross section is dependent upon the single variable r = 
Q 2 /Qs(x), as 

cr^*p(x,Q ) = <x-y*p(r). (4) 


3. Linear behavior for the evolution of the 
HFSFs 


The heavy flavor structure functions (HFSFs) are de¬ 
scribed as Mellin convolutions between the gluon distri¬ 
bution G(x,g 2 ) and the Wilson coefficients C k g {x, 0?) 
as 

Fk(x, Q 2 ) = Ctjx, %)®G(x,qi 2 ), (k = 2&L, i = c&b){ 5) 

I* 

here the Mellin convolutions is given by 

fi{x)^f 2 {x)= f — /i(a/)/ 2 ( —)- (6) 

J x V V 


The Q 2 evolution equation for the HFSFs (Eq.5) is ex¬ 
pressed in terms of these structure functions, as we have 


dFj:(x,Q 2 ) 

d In Q 2 


E 


dlna s 


n —1 

dlnC, 


d In Q 


^ n \x,Q 2 ) + P gg ®Fl{x,Q 2 ) 


+ ■ 


k,g 


dlnQ 2 


®F k {x, Q 2 ), 


( 7 ) 


in which the corresponding physical evolution kernel P gg 
can be derived from DGLAP evolution equation at small- 
x as it follows: 

dG dlnQ 2 ^ = 


This kernel is corresponding to the massless Wilson coef¬ 
ficients in leading order (LO) up to next-to-next-to lead¬ 
ing order (NNLO) [42-44], and also the heavy contribu¬ 
tions are in leading order and next-to-leading order by 
using massive Wilson coefficients in the asymptotic re¬ 
gion Q 2 » m 2 . 

According to the Regge pole approach, the distribution 
functions can be controlled by Pomeron exchange at 
small x, since these behaviors are correspondent to the 
BFKL (Balitskii-Fadin-Kuraev-Lipatov) Pomeron[45-48] 
ideas as extended by adding a hard Pomeron which de¬ 
scribes the small-a: HERA data up to Q 2 of a few hun¬ 
dred GeV 2 values. The small-x asymptotic behavior for 
the gluon and heavy flavors can be exploited to the evolu¬ 
tion equations of the HFSFs. Therefore, linear evolution 
of the HFSFs at small-a: can be found as 


F*(a:,Q 2 )Uo = E F fe (n) (^Qo)4 


( 9 ) 
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where 

ri,(n) 

1 k 


exp( / d In Q ( 
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( 10 ) 


Now, we have a compact small-a: formula for the heavy 

pi 

flavor ratio R l = which greatly simplifies the extrac- 
tion of F£ from measurements of reduced heavy cross 
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sections a‘ r : 

K(x,Q 2 ) = F^x,Q 2 )[ 1 - i+ y~_ y)2 R\x,Q 2 )\. (11) 


Therefore, we found a small-a; formula for the ratio R l 
by the following form 


R\x,Q 2 ) = 


FIM 2 ) _ Fi(x,Ql) 


Fl(x,Q 2 ) 


Fi(x,Q 2 ) 


exp( / din Q 2 ( 


d[\nCl g In Cl g ] 
d In Q 2 


( 12 ) 


<g>x A )). 


Here % = where Xo(— 0.01) is the boundary condition 
that the gluon distribution joints smoothly the unshad¬ 
owed region and R is the radius of the proton. The first 
term is the linear evolution (Eq.7) and the second term 
is due to the 2-gluon density. Therefore the shadowing 
corrections to the HFSFs can be defined by Eq.14. To 
obtain a differential equation for shadowing corrections 
to the HFSFs, we write out the sum and the coefficient 
function explicitly. Consequently we find an inhomoge¬ 
neous first-order differential equation which determines 
HFSFs shadowing corrections in terms of shadowed glu¬ 
ons. Eq. (14) can be rewritten in the following form: 


4. Nonlinear behavior for the evolution of the 
HFSFs 

As mentioned above the HFSFs linear evaluation equa¬ 
tion is based on a hard-Pomeron behavior for the gluon 
and heavy structure functions [49-56]. The gluon density 
increases with decreasing x and we must reach the region 
in which gluon-gluon interactions confine the growth im¬ 
plied by this behavior concerning the gluon distribution, 
G ~ x _A , as this behavior will violate unitarity when 
x— >0. Thus, we discuss absorption effects which tame 
the violation of unitarity. At sufficient small-x values two 
gluons in different cascades may interact and the gluon 
ladders fusion are generally important. Therefore the 
gluon density is decreasing and shadowing contributions 
can no longer be neglected [21-22,57-60]. Shadowing cor¬ 
rections, which take into account the fusion of t-channel 
gluons, modify the linear DGLAP equation for the gluon 
distribution by adding a negative term proportional to 
quadratic in g(x,Q 2 ). This picture allows us to write 
the GLR-MQ equation for the gluon distribution func¬ 
tion behavior at small-x symbolically as: 


dxg{x , Q 2 ) 
91n Q 2 


dxg{x,Q 2 ) 
dlnQ 2 


| DGLAP 


16f? 2 Q 2 


— [yg{y,Q 2 )} 2 - 
y 


(13) 


The nonlinear shadowing term, oc — [g] 2 , arises from per¬ 
turbative QCD diagrams which couple four gluons into 
two gluons so that two gluon ladders recombine into a 
single one. The minus sign occurs because the scatter¬ 
ing amplitude corresponding to a gluon ladder is pre¬ 
dominantly imaginary . Thus the equation (13) becomes 
nonlinear in xg. We neglect the quark-gluon emission di¬ 
agrams due to their little importance on the rich gluon 
in small-x region and we work under an approximation 
of neglecting contribution from the high twist gluon dis¬ 
tribution Ght{x, Q 2 ) [24]. 

In what follows it is convenient to use directly the re¬ 
duced gluon distribution function behavior according to 
the Eq.13 as the Q 2 evolution of the HFSFs modified by 


dFj:(x,Q 2 ) 
d In Q 2 


dFl(x,Q 2 ) 
9 In Q 2 


[Eq- 7] 
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where 


Sl(x,Q 2 )= 2e 
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and 


g(x,Q 2 ) = (J2 
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dlna s 

71 - 

1 d In Q 2 


Pgg®X X 


dlnC, 


k,g 


d\nQ 2 


<8>x x ). 


(17) 


The general solution of Eq. (15) where tames the behav¬ 
ior of the HFSFs at small-x has the following form: 


F z k (x,Q 2 ) 



Therefore, shadowing corrections to the HFSFs modify 
the heavy reduced cross section and also the ratio of the 
HFSFs, as we will have: 


KM 2 ) 


Fl(x,Q 2 )[Eq. 18][1- 1 + ( °_ y)2 
x&(x,Q 2 )[Eq.20]], (19) 


where 
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R l (x,Q 2 ) 


Fj(x,Q 2 ) 
Fi(x,Q 2 ) 


F[(X, Q 2 0 ) - /£ 2 Q 2 )ef ~VL(x,Q 2 )d\nQ 2 dln Q2 

r- -Fl -1 

Fi(x, Q 2 ) - f ®2 S z 2 (x, Q 2 ) e /-r /2 (x,Q 2 )dl n Q 2 d l n g 2 


x exp( / (ilnQ 2 ( 


d[lnCi jg - InC'l g] 

d In Q 2 


®x A )). 


( 20 ) 


Equations 18-20 satisfy the requirements expected for 
shadowed distributions of the heavy quarks. These 
equations reduced to the unshadowed distributions 
(Eqs.9-12) when shadowing is negligible; that is, when 
S l k (x, Q 2 )—>0 which implies G so t—>■oo. Finally the 
DGLAP+GLRMQ evolution joins smoothly into the 
DGLAP evolution when x—>xq- 


5. Geometrical scaling of the HFSFs 

In the saturation model [61], the dipole cross section is 
bounded by an energy independent value do (Eq.3) which 
imposes the unitarity condition, cr q q<ao with respect to 
the free parameters in the model[38-41]. These free pa¬ 
rameters can be extracted from data within some specific 
models of DIS. In the Golec-Biernat-Wiistlioff (GBW) 
model we see, a 0 = 23 mb , A = 0.3, Qo = 1 GeV/c and 
xq = 3xlCU 4 . The dipole cross section for a small qq 
dipole is related to the gluon density at scale g 2 as: 

<?qq = ^-r 2 a s xg(x,n 2 ). ( 21 ) 

For small r<^2Ro(x) in the saturation model (where 
Ro(x) is the saturation scale at small-x), the gluon den¬ 
sity found [38-41] by the following form: 


xg(x,g 2 ) 


3 cr 0 
47r 2 a s Rq(x) ’ 


( 22 ) 


where at the geometric scale for the boundary Q 2 = 
Q 2 (x) we have: 


^ x ff( x > Q 2 = Ql( x )) = r o x \ (23) 

with r 0 = -^ct 0 Xq. Also R.S.Thorne [62] used the re¬ 
lationship between the dipole cross section and the un¬ 
integrated gluon distribution and showed that the gluon 
distribution at fixed coupling has this behavior: 

xg(x,Q 2 ) = + (x 0 /x)* 

x(l - e - * 32 ^/® 0 ^)). (24) 


when a pair of heavy quarks HH (cc or bb) are produced 
in the final state. Therefore the gluon distribution can 
be evaluated for heavy quarks according to Eq.24. This 
is consistent with a general-mass-variable flavor number 
schemes (GM-VFNS) [63]. In this case, the parame¬ 
ters obtained from the best fit were (Jq = 29.12 mb , 
A = 0.277 and xq = 0.41 xlO -4 [38-41], Because the 

data for the bottom component are not suitable 

for scaling analysis( since they contain two few points 
[11-15]), therefore we used the charm parameters in our 
determination for bottom structure functions. 

Therefore the dipole picture at the geometric scaling is 
suitable for HFSFs analysis . If the characteristic size of 
the qndH ~Pair is much smaller than the saturation ra¬ 
dius, with decreasing x, it can show that the behavior of 
the HFSFs at the scale < g 2 > is given by 

Fk(x,Q 2 ) = r 0 Q 2 s (x)[Cl g (x, %)m x ]- (26) 

d 

We summarized our results in sections 3-5 at Figs.2-4. 
In these figures the results of calculation for charm 
and bottom structure functions, longitudinal structure 
functions and reduced cross sections are shown at 
Q 2 = 12 and 60 GeV 2 respectively. We determined the 
HFSFs and reduced cross sections in linear behavior 
(unshadowed) when shadowing is neglected. It was 
also made clear in nonlinear behavior (shadowed) at 
R = 5 GeV _1 where the gluons are spread throughout 
the entire proton and at R = 2 GeV _1 where the 
gluons are concentrated in hot-spots within the proton. 
We observe that these behaviors for the heavy quarks 
functions are tamed at small-ax The reduced cross 
sections are determined at the average inelasticity < y > 
in according with the HI data. We compared our results 
with GJR parameterization [64], the ZEUS and HI data 
[11-15]. We can observe that there are a well agreement 
between our unshadowed, shadowed and saturation 
results with the experimental data accompanied by total 
errors. 


In order to be able to study the formal heavy flavor pro¬ 
duction limit, the Bjorken variable x = Xb was modified 
into the one as follows: 

A_ rn ^ 

X = XB (1 + -Jf), (25) 


6. Heavy flavor exponents behavior 

In the double asymptotic limit, the behavior of the 
gluon distribution is expected to rise approximately as a 
power of x towards small-ax Because, the behavior of the 
HFSFs and heavy flavor reduced cross sections are depen- 
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dent on the gluon distribution behavior directly. There¬ 
fore we consider the power-like behavior of the HFSFs 
at small-a: as ^ where /= = F 2 ,Fl and er r . 

The logarithmic ^-derivative of the heavy flavor functions 
concerning ln(-) can be defined [65] as: 


H\{x,Q 2 ) 


din f{(x,Q 2 ) 
<91n(l/x) 


(27) 


Here we would like to consider the relation between 
the effective intercept A^ for heavy flavor functions and 
logarithmic-a; derivatives of the heavy flavor functions 
where 


1 dXUx, Q 2 




(28) 


For unshadowed heavy flavor functions (Figs.2-4) we used 
the Pomeron intercept in our determinations, therefore 
the effective intercept and (E-slop strictly coincide with 
the hard Pomeron intercept as we considered this behav¬ 
ior for the gluon distribution at small-a;. Therefore in 
this case, we will come into the below: 


H\(Q 2 ) = \\{Q 2 ) = 0.44. ( i = c,b ) (29) 


Figure 5 shows the derivative as a function of Q 2 for 
x = 0.001 value. The derivative is observed to have a 
logarithmically approximate rise with Q 2 for F%(x,Q 2 ). 
The ln(<3 2 ) dependence of F 2 (x,Q 2 ) is observed to be 
non-linear. It can be well described by a quadratic ex¬ 
pression, since for each of heavy quarks these derivatives 
can be described by the function b(x) + 2c(x) ln(Q 2 ). The 
shape of these derivatives reflect the behavior of the gluon 
distribution in the associated kinematic range. 

This scaling violation shows the transition from soft to 
hard dynamics. We give this scaling violation to the func¬ 
tion r/(Q 2 ) then we consider the heavy functions behavior 
as a fixed power of x: 

fl(x, Q 2 ) = v\{Q 2 )x~ X ■ (33) 

It would be implied that the Mellin transform f£(j,Q 2 ) 
would have a pole at j = 1 + A, where this singularity 
referred to as the hard-Pomeron singularity [49-56]. 
We conclude that the power A’s in Eqs.29 and 30 are 
not dependent on Q 2 and theses powers have got fixed 
values. In our determinations (Fig.5), each of theses 
coefficient functions r/((Q 2 ) vanished in a way that 
Q 2 — >0. In Fig .6 we compared our results for ^(Q 2 ) and 
rj 2 {Q 2 ) with DL model [49-56] and color dipole model 
(CDM) [1-4]. These results are comparable with others. 


When we consider the saturation effects [25-26] on the 
charm functions, the charm exponents obtained will be 
follows: 


HHQ 2 ) = A |(Q 2 ) = 0.277. (30) 

But this value is not coincident with the bottom func¬ 
tions, since we do not have enough data for bottom com¬ 
ponent in scaling analysis. One can see that exponents 
obtained for heavy flavor functions at the shadowed re¬ 
gion are dependent on the saturation scale and R pa¬ 
rameter. For the nonlinear evaluation equation we can 
observe that: 


H\{x, <3 2 )^A|(:e, Q 2 ). (31) 

For heavy functions at R = 5 and 2 GeU -1 , there is 
not a linear function when we fit it to all data because 
inequality (31) is correct. For R = 5 GeV~ 1 there is 
a second-order function and for R = 2 GeU -1 there 
is a three-order function. So we conclude that at non¬ 
linear evolution equations, the heavy exponents are de¬ 
pendent on x and R values, for the shadowed heavy 
flavors are dependent on the values x, Q 2 and R , i.e., 
(A { ’ =\{x,Q z 1 R)). 

In addition to the (E-slope, we want to consider the loga¬ 
rithmic derivative of the heavy functions with respect to 
InQ 2 , as the Q 2 -slopes are defined by: 


dfl{x,Q 2 ) 
dln(Q 2 ) ' 


( 32 ) 


7. Conclusions 

We have studied several aspects of the heavy flavors 
functions behavior at small-a;. To do this, we assumed 
that the heavy parton density obeys power-laws having 
effective power. Therefore we have applied the hard 
Pomeron and saturation methods in obtaining the heavy 
flavors functions considering HI and ZEUS data. We 
have shown that geometrical scaling in DIS works well up 
to the heavy flavors functions with a constant exponent. 
The merit of this exponent is mainly due to its relation 
with the gluon density at small- x. Results obtained 
suggest that geometrical scaling in heavy production is 
basically the same as for the inclusive DIS. The value of 
shadowed exponents are different from the one obtained 
for unshadowed corrections and geometrical scaling in 
heavy functions. This difference is due to the shadowing 
corrections on the heavy functions. Indeed the value of 
the shadowed correction on the exponents depends on 
how exactly the gluon ladders couple to the proton or 
on how the gluons are distributed within the proton, 
especially in hot-spot point. As a result, the coefficient 
functions ( ri(Q 2 )) are dependent on the Q 2 scale and the 
effective powers are constant. We obtained our results 
for the r](Q 2 ) and compared them with DL and CD 
models. 
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FIG. 1: Boson-gluon-fusion (BGF) graph. 
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FIG. 2: Heavy flavor structure functions F^{x, Q 2 ) as a func¬ 
tion of the Bjorken variable x at Q 2 values 12 and 60 GeV 2 
in comparison with the experimental data from ZEUS and 
HI Collaborations [11-15] accompanied by total errors and 
GJR parameterization [64]. The unshadowed and shadowed 
corrections and also the saturation model concerning the ge¬ 
ometrical scaling are shown in this figure. 
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FIG. 3: Heavy flavor longitudinal structure functions 
Fl(x, Q 2 ) as a function of the Bjorken variable x at Q 2 values 
12 and 60 GeV 2 . The unshadowed and shadowed corrections 
and also the saturation model concerning the geometrical scal¬ 
ing are shown in this figure. 
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FIG. 4: Heavy flavor reduced cross section a?.(x, Q 2 ) as a func¬ 
tion of the Bjorken variable x at Q 2 values 12 and 60 GeV 2 in 
comparison with the experimental data from ZEUS and HI 
Collaborations [11-15] accompanied by total errors at inelas¬ 
ticity < y >. The unshadowed and shadowed corrections and 
also the saturation model concerning the geometrical scaling 
are shown in this figure. 
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FIG. 5: Q 2 -slope heavy flavor structure functions 

(dF£/d\nQ 2 ) as function of Q 2 at x = 0.001. Our results 
shown for unshadowed and shadowed corrections and also for 
saturation model with respect to the geometrical scaling in 
Q 2 -slope heavy flavor structure functions. 
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